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Abstract. The observational link between Supermassive Black Holes (SMBH) 
and galaxies at low redshift seems to be very tight, and statistically the global 
evolution of star formation activity and BH accretion activity also seem to trace 
each other closely. However, pinning down the co-evolution of galaxies and BH 
on an object-by-object basis remains elusive. I present results from new models 
for the joint evolution of galaxies, SMBH, and AGN, which may be able to help 
resolve some of the observational puzzles. A unique aspect of these models is 
our treatment of self-regulated BH growth based on hydrodynamic simulations 
of galaxy-galaxy mergers. Although these models do quite well at reproducing 
the observed evolution of galaxies, they do not reproduce the observed history of 
BH accretion, predicting too much early accretion and not enough at late times. 
I suggest two possible resolutions to this problem. 



1. Introduction 

There is strong evidence that the present-day properties of Supermassive Black 
Holes (SMBH) and their host galaxies are tightly linked — for example we 
observe a tight relations hip between BH mass and the mass of the host spheroid 
in nearby galaxies (e.g. iHaring &: Rixl [20041 ) . However, we do not yet know 
whether this BH-galaxy mass relationship evolved over cosmic time, or how 
galaxies and SMBH arrived onto it. Perhaps most importantly, we do not yet 
understand the physical origin of this relationship. 

Mergers have often been proposed as a mechanism that can drive gas 
into the nucle i of ga laxies and fuel both starbursts and black hole accretion. 
Hopkins et al.l ( 20051 ) used a large suite of hydrodynamic simulations of galaxy 



mergers to characterize the e pisode of black hole growth and the AGN accretion 
"lightcurves" in these events. iHopkins et al.l (p006 ) then used this model to test 
whether observed galaxy merger rates, the build-up of red/spheroidal galaxies, 
and AGN luminosity functions are consistent with the picture in which mergers 
both trigger AGN activity and cause morphological and color transformation. 
They concluded that there is excellent statistical agreement between all of these 
quantities. 

However, establishing a direct link between mergers and AGN activity has 
proven controversial. Low luminosity AGN, which are more common and can 
be studied in more detail at relatively low redshift (z < 1), do not show strong 
evidence f or being predominantly in morp hologically disturbed hosts or close 
pairs (e.g. iLi et al. l2006l : iPierce et al.ll2007l ). Higher luminosity quasars seem to 



be more often (th ough not always) associated with disturbed hosts or to hav e 
close companions (iBahcall et al.lll997l : iBennert et al.l l2008l : iLetawe et al.l 120081 ') . 
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but the samples with high resolution imaging are small. It may be the case that 
AGN are most easily identified in the very late stages of mergers, when the two 
nuclei have completely coalesced and most signs of morphological disturbance 
have become invisible. Especially at high redshift, at the typically available 
image sensitivity, these very late stage merger remnants might well appear to 
be "normal" spheroidal- type galaxies. 

Of course it is also possible, even likely, that there are multiple modes of 
BH growth. There are ample suggestions in the literature that AGN activity 
could also be fed by b a rs or b y stochastic accretion of cold gas in galactic nuclei. 
Hopkins Sz Hernquist ( 20061 ) presented a model for fueling of BH in isolated 



disks by stochastic accretion, and argued that while this fueling mode may be 
important in lower luminosity objects and at late times, most of the growth of 
today's massive BH occurred in the merger-driven mode at high redshift. 

Another sign of co-evolution is that the global histories of star formation 
and BH accretion seem to trace each other remarkably closely over the whole 
range of lookback times over which these quantities have robust observational 
estimates. The global star formation rate density is almost exactly a factor of 
2000 times the BH accretion rate density, from redshift z ~ 5 to the present. 
Moreover, even when divided by mass or accretion rate, "matched" populations 
of galaxies and black holes at z < 1 also trace one another's activity (Zheng et 
al. in prep). However, we know that at least at z < 1, the bulk of the star 
formation activi ty is associated with isolated disks, not mergers or spheroids 
( Bell et alll2005l ). Therefore, it is also strange that these two kinds of activity 
trace each other, but do not seem to be occuring in the same types of objects 
(Zheng et al. in prep). 

At higher redshi ft, there may be a larger mis match between between the 
two kinds of activity. iFaucher-Giguere et al. ( 20081 ) obtained constraints on the 
global star formation density at 2 < z < 4.2 from the Lyman-a forest opacity 
in QSO spectra. They found that the hydrogen photoionization rate, and hence 
the star formation rate density, was remarkably flat over this redshift range, in 
contrast to the sharply peaked QSO luminosity density, which falls off sharply 
at z > 2. This may indicate that there is a time delay between SF activity and 
QSO activity, which is naturally explained in the merger picture (any galaxy 
with cold gas can form stars, but galaxies have to "wait" for a major merger 
before significant accretion onto the BH can occur). 

One way to test this picture, in which mergers are responsible for triggering 
both AGN activity and morphological and spectrophotometric transformation 
of galaxies, is to build cosmological models that treat the growth of galaxies, 
black holes, and AGN self-consistently. In addition, many astronomers now 
believe that the energy released by accreting black holes may play a crucial role 
in regulating galaxy formation. Here we describe one such semi-analytic model 
for the joint formation of galaxies, black holes, and AGN, and present some 
predictions from these models. 



2. Unified Models for Galaxies, Black Holes, and AGN 



The foundations of our model, which pertain to the growth of structure in 
the dark matter component, and the formation of galaxies, are described in 
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Briefly, the models are implemented within dark matter halo "merger trees" 
and include approximate treatments for atomic cooling of gas, photoionization, 
the formation of angular momentum supported disks (including estimates of 
their sizes and circular velocities), star formation according to an empirical 
"Kennicutt" Law, supernova feedback, and chemical evolution. We have re- 
cently implemented new machinery within this framework to treat black hole 
gr owth and the ass ociated AGN feedback. Our new models are fully described 
in ISomerville et"al] (|2008 . SOS). Here we give a very brief synopsis of the most 
important new ingredients. The models presented here are for a "Concordance" 
ACDM cosmology with Qrn = 0.3, 17a = 0.7, Hq = 70 km/s/Mpc, as = 0.9, 
n.o = 1. 



Galaxy-galaxy Mergers We assume that galaxy mergers may induce a burst of 
star formation and destroy any pre-existing disk component, depending on the 
mass ratio of the smaller to the larger galaxy, fi. We parameterize the strength 
and timescale of these bursts according to the results of a large suite of hy- 
drodynamic simulations of galaxy mergers ( Robertson et al. 2006b : ICox et ahl 
l2008l ). as described in SOS. Mer gers can also heat and thicken, or even destroy, 
a pre-existing disk component, driving galaxies towards morphologically earlier 
Hubble types. We assume that the fraction of the pre-existing stars that is trans- 
ferred from a disk to a spheroidal component is a strongly increasing function of 
the merger mass ratio ^, such that minor mergers with /j, < 0.2 have little effect, 
and major mergers with fj, > 0.25 leave behind a spheroid-dominated remnant. 



Bright Mode AGN and AGN-driven winds In our model, mergers also trigger 
the accretion of gas onto supermassive black holes in galactic nuclei. Each top- 
level halo in our merger trees is seeded with a black hole of mass Mgeed — 
IOOMq. Following a merger, the black hole is allowed to grow at the Eddington 
rate until the BH reaches a critical mass, where the radiative energy being 
emitted by the AGN becomes sufficient to halt further accretion. This self- 
regulated treatment of black hole g rowth is based on hydrodynamic simulations 
including BH growth and feedback ( Springel et al.|[2005bl : iDi Matteo et al.ll2005l : 



IS , , 

Hopkins et al.ll2007al ). and is described in more detail in SOS. Energy radiated by 



black holes during this "bright", quasar-like mode can a lso drive galact ic-scale 
winds, clearing cold gas from the post-merger remnants (jSpringel et al.ll2005a ). 
Our model for momentum-driven AGN winds is described in SOS. 



Radio Mode Feedback In addition to the rapid growth of BH in the merger- 
fueled, radiatively efficient "bright mode" described above, we assume that BH 
also experience a low-Eddington-ratio, radiatively inefficient mode of growth 
associated with efficient production of radio jets that can heat gas in a quasi- 
hydrostatic hot halo. The accretion rate in this phase is m odelled assuming: 
Bondi accretion using the isothermal cooling flow solution of iNulsen Sz FabianI 
(|2000l 1. We then assume that the energy that effectively couples to and heats 
the hot gas is given by Lheat = «^heat^"^radioC^, where mj-adio is the accretion 
rate onto the BH, r/ = 0.1 is the assumed conversion efficiency of rest-mass into 
energy, and Kheat is a free parameter of order unity. 
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Figure 1. The predicted relationship between bulge mass and black hole 
mass (grey shading indicates the conditional probability -P(mbhl'Tibuige); light 
blue solid and dashed line shows the median an d 16th an d 84th percentiles) 
compared with the observed relation from Har ing fc Rixl (1200 4, green lines). 
Symb ols show the measurements for individual galaxies from iHaring fc Rixl 
(l200l . 



3. The Global Evolution of Galaxies and Black Holes 

We normalize the free parameters in our models to reproduce a key set of ob- 
servational quantities at z ~ 0, such as the stellar mass function, the cold gas 
fraction in disk galaxies, and the stellar metallicity (SOS). We show our model 
predictions and compare with observations for a broad range of low redshift 
galaxy properties in SOS. 

The relationship between galaxy mass and BH mass is clearly a key result 
that our model should reproduce. In our model, this relationship is set by the 
depth of the potential well of the galaxy at the time when the BH forms, which 
in turn is determined by the gas fraction of the progenitor galaxies of the last 
merger. More gas-rich progenitors suffer more dissipation when they merge, 
and produce more compact remnants with deeper potential wells. A deeper 
potential well requires more energy, and therefore a more massive BH in order 
to halt further accretion and growth. We see from Fig. [1] that our fiducial model 
reproduces the observed slope and scatter of the MBH--^sph (black hole mass 
vs. spheroid mass) relationship. Our model also predicts that the MeH^-^sph 
relation should evolve with time. Because more gas-rich merger progenitors 
produce remnants with larger black holes, and the galaxies in our models were 
significantly more gas rich in the past, galaxies have larger BH for their spheroid 
mass than they do at the present day. This leads to a relatively mild amount 
of evolution, of a factor of less th an two since z ~ 1 and a factor of about four 
since z ~ 3 ( Hopkins et al.l [2"007al ) . 



Fig. [2] shows the global star formation rate density of all galaxies predicted 
by our models. We show both the results of our "fiducial" model, as well as 
a model in which galaxy formation has been suppressed in low-mass halos (we 
simply do not allow gas to cool in halos with mass less than 10"^"^ Mq). We will 
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Figure 2. Left: Star formation rate density as a function of redshift. The 
upper set of thicker lines shows the total SFR in the models, and the lower 
set of thin lines shows the SFR due to bursts. Symbols and solid lines show 
a compilation of observational results, converted to a Chabrier IMF (see SOS 
for details). Right: The integrated global stellar mass density as a function of 
redshift. Symbols and the solid line show observational estimates (see SOS). 
In both panels, dashed (purple) curves show our "fiducial" ACDM model, and 
dot-dashed (red) curves show a model with reduced galaxy formation in small 
mass halos (see text). 



refer to the latter model as the "low" model for brevity. The results of the low 
model are similar to those from models with reduced small scale power, e.g. with 
lower (Tg. We have kept the same values for the free parameters in both models. 
Both models fit the data fairly well from < z < 1, and are about 0.1 5-0.2 
dex lower than the observational compilation of Hopkins Sz Beacom ( 20061 ) from 
1 < z < 3. The fiducial model predicts much more star formation at very high 
redshift {z > 3) than the low model, because much of the star formation at these 
redshifts is taking place in low-mass halos. 

In the right panel of Fig. [2] W6 show tliG complGmGntajry qua^ntity /^star? 

the 

integrated cosmic stellar mass density. The fiducial model predicts a significantly 
earlier assembly of stars in galaxies than observations of high redshift galaxies 
indicate. However, the low model produces very good agreement with the stellar 
mass density as a function of redshift. We note that this tension in the model 
results is connected with a possible inconsistency between the two observational 
data sets (star forma tion rates and stellar masses) that has been noted recently 
in several papers (e.g. lHopkins Beacom|[2006l : iFardal et al]l2007l : IWilkins et al 



20081 : Pavel [200^ ). One possible resolution of this tension can be obtained if the 



stellar IMF has changed with time, and was more top-heavy at high redshift. 



4. Evolution of "Bright Mode" AGN 



Fig. [3] shows the integrated bolometric luminosity density of "bright mode" BH 
accretion in our mo dels, compared with o bservational estimates from quasar lu- 



minosity functions ( Hopkins et al. 2007bl ). We see that in our fiducial model. 



QSO activity peaks at too high a redshift, and falls off too rapidly at low red- 
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Figure 3. The bolometric luminosity density of "bright mode" AGN 
(quasars) as a function of reds hift. Symbols w ith error bars show the obser- 
vational estimates from Hopki ns et alJ ()2007bl ) . Lines show model predictions 
from: our fiducial model (solid), our "low" model (short-dashed), a model in 
w hich black holes are e jected after every merger (long dashed), and the model 
of iCroton et all (p006l) . 
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shift. The "low" model predicts a rise and fall in QSO activity that is in some- 
what better agreement with the data, but still does not have quite the right 
shape. For compari son, we plot the bright mode accretion rate in the models of 
(jCroton et al.ll2006l ^. which are similar to our fiducial models in many respects, 
but treat BH growth and AGN feeding differently. We see that this model also 
fails to reproduce the observations: QSO activity peaks too early, and does not 
fall off rapidly enough at low redshift. 

This result is surprising, because iHopkins et aD (|2008 b"l found good agree- 
ment with these same observations, using the same model for BH growth that 
we have implemented here. We checked that the merger mass function predicted 
by our model agrees with that in the Hopkins et al. models. The gas fractions 
of our merger progenitors also agree well. E ventually we realized t hat an im- 
portant difference in the assumptions made bv lHopkins et al.l (|2008bl 1 and those 
contained in the models presented here was that, in the Hopkins et al. model, 
the black holes were allowed to grow from the seed mass to the final mass in 
each merger. This effectively neglected the pre-existing BH that grew in earlier 
mergers. Put another way, in our models, by 2 ~ 1-2, there is already a sig- 
nificant amount of mass in spheroids in place, each containing a massive black 
hole. If the pre-existing BH is already above the critical mass discussed above, 
then no AGN activity takes place. Thus, in our models, many late mergers do 
not trigger any AGN activity. 

As a check, we re-ran our model, assuming that the black holes are ejected 
from the galaxy at the end of every merger. The bright mode accretion history 
now agrees much better with the observations at least in shape, although the 
normalization is a bit high. Although this model may seem artificial, it is possible 
that in some cases black holes are actually ejected by the gravitational rocket 
mec hanism, which m ay impart kicks of several hundreds to 1000 km/s to merging 
BH ( Volonteril 12007 . and references therein). However, this ejection of massive 
BH can clearly have a problematic impact on the predicted MBH~-^sph relation, 
as shown in Fig. reffig:mbh. However, the model shown here is extreme in that 
BH are ejected from every galaxy after every merger. It is possible that a more 
physical implementation of the gravitational rocket mechanism could improve 
the BH accretion history without leaving behind too many massive spheroids 
with no BH, or with very low mass BH, which are not observed. 

Another possible resolution to this problem could perhaps be obtained if 
spheroids are not formed as efficiently in high redshift major mergers as we have 
assumed here. Because in our models, the properties of the spheroid regulate 
BH formation, if spheroids form later, then BH accretion will occur later as well. 
There are indications from a detailed analysis of the hydrodynamic merger simu- 
lations th at gas-rich mergers may often produce disk-like , rather than spheroidal, 
remnants ( Robertson et al.ll2006al : [Hopkins et al.ll2008al ) . Because our disks tend 
to be more gas rich at high redshift, this effect would indeed shift spheroid for- 
mation to later times. We are currently working on implementing this effect in 
our models. 
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5. Summary and Conclusions 

We have presented some first predictions from new models for the joint evolution 
of galaxies, black holes, and AGN. Our models differ from other semi-analytic 
models in the literature in that we have implemented self-regulated BH growth 
based on the results of high resolution hydrodynamic simulations of galaxy merg- 
ers. 

We showed that our fiducial model does a fairly good job of reproducing the 
global star formation history, but overproduces the stellar mass in place at high 
redshift. A "low" model with cooling suppressed in low-mass halos, or a model 
with reduced small-scale power, produces good agreement with the stellar mass 
density history but does not agree as well with the star formation history. This 
reflects an internal tension between these two data sets, which may be due to 
systematic errors arising from a time- varying stellar Initial Mass Function. 

We found that our model, which does fairly well at reproducing galaxy 
observations at both low and high redshift, fails to reproduce the form of the 
BH accretion history as traced by the bolometric quasar luminosity density. In 
our models, QSO activity peak s at too high a redshift, and falls off too rapidly 
at low redshift. The models of ICroton et al. ( 20061 ) also fail to reproduce this 



behavior. We suggested two possible resolutions to this problem: 1) Black Holes 
are (at least sometimes) ejected from their host galaxy following a merger or 
2) major mergers at high redshift are not as efficient at producing spheroidal 
remnants as we have assumed. 
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